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Abstract. A field-theory approach is used to investigate the "spin-glass effects" on 
the critical behaviour of systems with weak temperature-like quenched disorder. The 
renormalization group (RG) analysis of the effective Hamiltonian of a model with 
replica symmetry breaking (RSB) potentials of a general type is carried out in the 
two-loop approximation. The fixed-point (FP) stability, recently found within the 
one-step RSB RG treatment, is further explored in terms of replicon eigenvalues. We 
find that the traditional FPs, which are usually considered to describe the disorder- 
induced universal critical behaviour, remain stable when the continuous RSB modes 
are taken into account. 



PACS numbers: 64.60.Ak, 64.60.-i, 64.60.Fr, 75.50.Lk 
1. Introduction 

The spin-glass effects on the critical behaviour of weakly disordered systems have 
received intense attention during the last several years [1-7]. It is well-known that 
the replica method is one of the very few general analytical methods to investigate 
disordered systems |5]. It enables one to average over disorder and employ standard 
field theory techniques, such as the RG analysis. Most of the results obtained using 
the RG methods implicitly assume replica symmetry (RS). However, as is well-known 
from the large body of works on glass-like systems, the Parisi RSB scheme is required to 
obtain the correct low temperature solutions of random models which are characterized 
by a macroscopic number of ground states [HI El- Recently, Le Doussal and Giamarchi pQ 
pointed out that although the RG treats fluctuations exactly, it has to be incorporated 
with RSB to escape a risk that it will miss the physics associated with the existence of 
a large number of metastable states in certain disordered systems. 

Dotsenko et al [2-4] studied the RSB effects on the critical behaviour of disordered 
p-component magnets to one-loop order using e expansion. They found that for p < 4 
the FPs, which are usually considered to describe the disorder-induced universal critical 
behaviour, are unstable with respect to the introduction of RSB potentials. The one-step 
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RSB ansatz 8j allowed them to find the new FPs which turned out to be stable within 
the one-step RSB subspace, but the further analysis showed that for a general type of 
RSB there exists no stable FPs and the RG flows lead to the strong coupling regime [3]. 
However, the numerous investigations of pure and disordered systems performed in the 
two-loop and higher orders of the approximation prove that the predictions made in the 
lowest order of the approximation, especially on the basis of the e expansion, can differ 
strongly from the real critical behaviour [TO]. Therefore, the results of the RSB effects 
investigation in [2-4] must be revised with the use of more accurate approximations. 

In |2j the one-step RSB effects on the critical behaviour have been reconsidered 
using the field-theory approach in the two-loop approximation. It was shown that the 
RS FPs are stable with respect to the introduction of the one-step RSB potentials. 
However, it is likely that the RSB in disordered systems should be of a general type 
rather than in the one-step RSB class Therefore, one cannot be sure that FP 
scenario obtained by means of the one-step RSB ansatz describes adequately the spin- 
glass effects on criticality. In the present work, we re-investigate the RSB effects on the 
critical behaviour extending the field-theoretic description given in |H] to the case of the 
RSB structure of a general type. 



2. Model and renormalization 



The starting point of the field-theory approach to the study of spin systems near 
critical point in the presence of weak quenched disorder is the Ginzburg-Landau- Wilson 
Hamiltonian 



H 



[ d d r Ul><Mr)] 2 + 1^1- *r(0] £tf (0 + (r)#(r)l (1) 

J {. i=i ' i=i i,j=i ) 

where <f>i(r) is the p-component order parameter. The quenched disorder is described 
by random fluctuations of the effective transition temperature 5r(r) whose probability 
distribution is taken to be symmetric and Gaussian 

1 

where u is the small parameter which is proportional to the concentration of defects; 
and A is the normalization constant. To carry out the averaging over disorder we use 
the standard replica trick. After replicating, we get the effective Hamiltonian 

p 71 P 71 P 71 



P[5t] = Aexp 



d d r(Sr(r)f 



(2) 



i=l a=l i=l a=l i,j=l a, 6=1 

which is a functional of n replications of the original order parameter with an additional 
vertex u in the replica symmetric matrix v a i = v8 a i, — u. The properties of the original 
disordered system are obtained in the replica number limit n — ► 0. 

As pointed out in j2j, near the critical point the system has a macroscopic 
number of the local minima solutions of the saddle-point equation corresponding to the 
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Hamiltonian (|T|). The traditional RG treatment of the replicated Hamiltonian (jHJ) cannot 
take into account the existence of these local minima so that the direct application of 
the RS RG scheme may be questioned. It was argued that spontaneous RSB can occur 
due to the interaction of the fluctuating fields with the local non-perturbative degrees of 
freedom coming from the multiple local minima solutions of the saddle-point equations. 
It was shown that the summation over these solutions in the replica partition function 
can provide the additional non-trivial RSB potential b tw>0 2 2 in which the matrix 
v a b has the Parisi RSB structure [Hj. 

We now realize the field-theoretic RG analysis of the effective Hamiltonian Q 
with the RSB potential in the two-loop approximation. We use the massive field 
theory scheme with renormalization of the one-particle irreducible (1PI) vertex functions 
Iq (A4, .., Ul\ <Zi, ••, ?iv; A*oj { v ab}) at non-zero mass and zero external momenta [TT] . 
The 1PI vertex functions are defined by 

5 (E^ + E^') r o L ' JV) ({ A; }'W^o,K b }) = / e^+<^) 

x (0(n) . . . 0M0 2 (ri) • • • 2 (r L )) 1PI d d n . . . d%d d f x . . . d d f L (4) 

where {q}, {k} are the sets of external momenta, the brackets (. . .) denote the averaging 
performed with the effective Hamiltonian (jHJ), and for the sake of simplicity the replica 
and vectors indexes have been omitted. The Feynman diagrams that contribute to the 
vertex functions involve the momentum integration and depend on the microscopic 
cut-off A . To study the critical domain which corresponds to the large cut-off limit 
(Ao — > 00) we have to renormalize the theory in order to absorb the divergences 
of diagrams in a change of parameters and to obtain meaningful expressions for the 
correlation functions. Introducing the renormalization factors for the fields Z^, Z^, the 
renormali2 
as follows, 



renormalized vertex functions T^' N ^ are expressed in terms of the bare vertex functions 



^ N \{k}; {<?}; /i 2 , {9a*}) = Z N J 2 ~ L Zy^ N \{k}- {q}- & {v ab }) (5) 

where fi is the renormalized mass and g a b is the renormalized matrix of couplings. To 
define the renormalization scheme completely, one has to impose the renormalization 
conditions for the renormalized vertex functions: 



k=0 



r ( R 2) {k,-k;Li 2 ,{g ab }) 
^rg' 2) (/c, -fc;/i 2 , {g ab }) 

rS'i ^i' fc 2, fa, h ; v 2 , {gab}) 

r^' 2) (/c, -k; g;/i 2 , {g ab }) 



l> 2 



k=0 



ki=0 

= 1. 



/i 4 d g a b 



(6) 



k=q=0 

The scaling behaviour in the critical domain is described by the homogeneous Callan- 
Symanzik equation for the vertex functions [TT] : 

+ ^ Pa' b' ({gab}) 



a 
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L l<f>({9ab})-L^^({g ab }) 



R 



({k};{q}- f i 2 ,{g ab }) = (7) 



with the coefficients given by 
dg a b 



Pab({gab}) 



lA{9oh}) 



dZ, 



{Vab}, m> 



<91n/i 



(8) 



{Vab},HO 



: d In u 

{VabhVO " 

Equations (jHJ) imply that the renormalized matrix g ab has the Parisi RSB structure, 
which is stipulated by the RSB structure of the bare matrix v ab . According to the 
technique of the Parisi RSB algebra, in the limit n — > the matrix g ab is parameterized 
in terms of its diagonal elements g and the off-diagonal function g(x) defined in the 
interval < x < 1: g a b —> (g,g(%))- The rules for operations with matrices g ab are 
detailed in |U|H|- in |S] the renormalization procedure ©-© was carried out with the 
assumption that all matrices g ab have the structure known as the one-step block-like 
RSB with the function g(x) having the form 

{On for < X < Xn 

(9) 
gi for x < x < 1 

where the undetermined step parameter 2 6 [0, 1] is related to the initial disorder 
distribution. In the present work, we renormalize the theory to two-loop order assuming 
that the RSB structure is of the general Parisi type. The renormalization conditions © 
can be expressed in terms of g and g(x) as follows 



[k } -k;^ ,g,g(x)} 



k=0 



r (0,2) 
1 R 

7^ r 2' 2) [k,-k;fi 2 ,g,g(x)] 



k=0 



=.(0,4) 
R 



[k u k 2 ,k 3} h;ij ,g,g(x)] 



(10) 



r i?' 4) i k u k 2 , h, k 4 ; g, g(x)] 



fe l =o 



T { ^\k,-k]q]n\g,g(x)} 



1. 



9\F) 



k=q=0 



Equations (fTUj) are simultaneous integral equations, which show that we deal with the 
functional RG. To obtain the renormalized couplings and the renormalization factors 
we can solve equations (jlU|) perturbatively order by order in g and g(x). The functional 
counterparts of equations (jHJ) read 



dg 



d In /i 



dg(x) 



'.ID 



S,«(*),W dln ^ 

where (v,v(x)) represents the bare matrix v ab . Equations (|TT|) allow us to calculate the 
(3 functions as functionals of the renormalized parameters g and g(x). In the two-loop 
approximation we obtain the following expressions for the (3 functions, 

£&0(aO]=7i+T2 + 73 ( 12fl ) 
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73 = ((8 / - 40 h + 20)p + 16f-176h + 88)g 3 + Ap(6h-2f-3)gg^ 
+ 8p(l-2h)g~ 3 ~ 

P\g,g(x)] = 7l + 72 + 73 (12b) 
7i = -f ( s )> 72 = -4fi- 2 (x) + (4 + 2p)gg(x) + 2pg(x) + pRi(x) 
73 = 16(1 - 2h)g 3 {x) - (8p(l - 4/i) + 48 - 96h)gg 2 (x) + ((8/ - 48/i + 28)p 
+ 16/ - 48/i + 24)^(x) + 16p(l - 2h)gg I g{x)+8p(l - 2h)gR x (x) 
- 8p(l - Ah)g^g 2 (x) + 8phg(x)R 1 (x) + 4p(4/i - 2/ - 3)^(x) 
-8p(l-2/i)i2a(x) 

1 X X 

where we have introduced g™ = J g n (y)dy, Ri(x) = J[g(x) —g(y)] 2 dy, R 2 {x) = f[g(x) — 

oo o 
g(y)][g 2 (x) — g 2 (y)]dy. The values of two-loop integrals are given by f(d = 3) = JL 

h(d = 3) = |, /(d = 2) = 0.11464, /i(d = 2) = 0.78129. Analogous to [2-5], we changed 

Qab —* gab/ J with J = f d d k/(k 2 + l) 2 and g a ^ — >■ —g a ^b in equations (12). The RS 

situation corresponds to the case g(x) = go in which equations (12) reduce immediately 

to equations obtained previously in |T2| for vertices v% = (p + 8)(g + go) and v<i = 8go. 

It is well known that the expansions of RG functions in powers of coupling constants 

are asymptotic series. A self-consistent way to extract the required physical information 

from the obtained expressions (12) requires application of special resummation methods: 

Borel resummation accompanied by certain additional procedures (see reference ^S] 

and citation therein). We employ the Pade-Borel resummation method extended to the 

functional case. The resummation procedure consists of several steps: (i) starting from 

the functional f[g, g(x)] in the form of a series in the auxiliary variable 8, 

f[g,g(x)) = Y,^ dm ^Y,^ dm ( 13 ) 

lit. 

m m 

where 7 m denotes the contribution of order g m , we construct its Borel transform (J13)) : 
(ii) the Borel transform is extrapolated by a rational Pade approximant [K/L\{9) that 
is defined as the ratio of two polynomials both in the variable 9 of degrees K and L such 
that the truncated Taylor expansion of the approximant is equal to that of the Borel 
transform of the functional /; (iii) the resummed function f s is then calculated as the 
inverse Borel transform of this approximant: 

oo 

fs[g,g(x)\ = J d6exp(-6)[K/L}(6). (14) 
o 

In the two-loop approximation (3 functions (12) have the form / = ^iQ + 72# 2 + 73# 3 . 
Using the [2/1] Pade approximant [2/1] (6) = ^6 + 37 2 # 2 /( 6 72 - 27 3 #) we obtain 

f\~ CM 37 ' 9 7 2 3 27 72 4 / 3t 2 \ . /3 72 \ 

fs[g, g{x)\ = 7i - —J + -77-3- exp Ei 15 

2 73 2 7 2 2 7 | V 7s / V 7s / 

where Ei(x) = V.V. dtexp(t)/t is the exponential integral. 
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3. Stability of FPs and replicon eigenvalues 

The nature of the critical behaviour is determined by the existence of a stable FP 
satisfying the following simultaneous equations 

P,\$*,9*{x)]=0, (3 s [g*,g*(x)}=0. (16) 

Taking the derivative over x, we obtain from the last equation 

^(3 s [g,g(x)} = X[g,g(x)]g'(x). (17) 

Equation (JT7jl implies that either the FP function g*(x) has the step-like structure, or 
it obeys the equation X\g*,g*(x)] = 0. The simplest assumption is that the one-step 
RSB (jHJ) occurs in the weakly disordered systems [2] so that the (3 functionals (12) 
reduce immediately to the three functions: Pi(g, go, gi), (i = 1,2,3) obtained in 0. 
The analysis of the (3 functions shows that there are three types of non-trivial FPs for 
different values of p — 1, 2, 3. Type I with g* ^ 0, g$ = g{ = corresponds to the RS FP 
of a pure system, type II with g* ^ 0, g^ = g{ ^ is a disorder-induced RS FP, and type 
III with g* 7^ 0, g£ = 0, g* ^ corresponds to the one-step RSB FP. The parameters 
9*,9o,9i, describing the FPs position, depend on the coordinate of the step x G [0, 1] 
and the obtained values of these parameters have been presented changing with the step 
Ax = 0.1 in tables 1-4 of paper 0. 

The type of critical behaviour of this disordered system for each value of p is 
determined by the stability of the corresponding FP. Within the one-step RSB parameter 
subspace, the requirement that the FP be stable reduces to the condition that the 
eigenvalues A» of the matrix 



B 



" " % 



Jj 

lie in the right complex half-plane. The values Aj calculated in [5, show that for three- 
dimensional (3D) and two-dimensional (2D) Ising models (p = 1) and 3D XY model 
(p = 2) the disorder-induced RS FPs are stable in the one-step RSB subspace, whereas 
the critical behaviour of 3D Heisenberg model (p = 3) is described by the RS FP of a 
pure system. 

As shown in |J] the stability of the FP in the full RSB space may strongly differ 
from the stability in the one-step RSB subspace, so it is quite interesting to investigate 
the stability properties of the calculated early FPs with respect to the RSB modes of a 
general type. The (3 functions (12) obtained in the present work give the complete RG 
description of system © within the scope of the two-loop approximation. In order to 
investigate the full stability properties of FP we have to solve the RG flow equations 
d 

J (19) 



Replicon modes and stability of critical behaviour 



7 



in the vicinity of the FP for // — > [11 . This task seems to be very complicated in 
the two-loop approximation. Recently, an alternative method, the so-called "replicon 
modes" approach, was suggested for this complex problem [§]. This approach allows 
for a very simple way to explore the stability properties with respect to the continuous 
RSB modes without solving equations ([!§]) directly. Indeed, the usual way to examine 
the stability of the FP of a model with several coupling constants is to linearize and 
diagonalize the RG flow equations about this FP. The replicon modes approach is the 
generalization of this method to the specific system (JHJ) with the infinite number of 
coupling constants. Following works E] we introduce the eigenmode function 

Q(x) = ±\g(x)-g*(x)] (20) 

which characterizes the behaviour of the first derivative of g (x) near the FP in the full 
RSB parameter space. In the case of a step-like RSB FP (RS FP, one-step or many-step 
RSB FP), as it follows from equation (|17|h the flow equation for the eigenmode function 
has the form 

d_ 

' d\i 

where we have introduced the replicon eigenvalues A rep (x) determining the scaling 
behaviour of the replicon modes Q(x). Thus, the change from g(x) to Q(x) diagonalizes 
the flow equation about the step-like RSB FP. If A rep (:c) is positive, Q^x) decreases 
under the iteration of the RG transformation and the function g(x) will approach the 
FP function g*(x). If A rep (a;) is negative, g(x) will deviate from its FP expression with a 
related instability along the continuous replicon mode directions in the parameter space. 
If A rcp (x) = 0, g(x) will keep the initial shape and this should correspond to a marginal 
stability. 

Using equations (fiTj|) - (|2*Tj) we derive A rep (x) = X[g* , g*(x)], where X[g,g(x)} is given 
by equation (JTTj). Substituting into the latter the resummed (3 function (|T5|) with 
the coefficients ji taken from equation (j!2fcj) we obtain the expression for the replicon 
eigenvalues function 



H—.Q^) = Kep(x)Qfi( X )i ( 21 ) 



A rep (x) = 



m=l 



1 <97m~ 


df s 


g'(x) dx 


{g*,g*(x)} d lrn 



(22) 

{ii=ii[g* ,g* (x)]} 

The set A = {A 1; A 2 , A 3 , A rop (x)} composed of the replicon eigenvalues and the 
eigenvalues of the stability matrix ()18j) yields a complete description of the FP stability 
properties in the full RSB parameter space [7j. We have calculated the replicon 
eigenvalues numerically for all FPs, which were found in the one-step RSB parameter 
subspace. Although the coefficient 73 in equation (|126|) depends explicitly on x, it turns 
out that in the case of the one-step RSB FP the replicon eigenvalues function (j2^j) 
depends on x through g*(x) only, and therefore, it also has the one-step structure 

f a££ for < x < x 
A rep (*) = . - ° (23) 

[ A; ep for x < x < 1. 
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Table 1. Replicon eigenvalues X, 



Type 


Xq 


d = 3 






d= 2 


p=l 


p = 2 


P = 3 


p=l 


T 

1 




a i f a o o f 

—0.169236 


—0.001673 


0.131538 


A AAn OA 

—0.095180 


TT 

11 




0.211760 


0.000003 


A A i < rni 

0.014501 


A A f A r A 

0.056959 


III 


0.0 


0.211760 


0.000003 


A A i < rni 

0.014501 


a a r o a r a 

0.056959 




0.1 


0.245366 


0.000092 


A AAAAAA 

0.009909 


A A/ 1 Z 1 1 AO 

0.066108 




0.2 


0.283705 


A A A A 1 A 1 

0.000191 


A AA rnin 

0.005010 


0.076207 




0.3 


0.326606 


0.000302 


-0.000202 


0.087391 




0.4 


0.375427 


0.000425 


-0.005714 


0.099813 




0.5 


0.430949 


0.000566 


-0.011485 


0.113656 




0.6 


0.494819 


0.000728 


-0.017414 


0.129128 




0.7 


0.568823 


0.000913 


-0.023284 


0.146478 




0.8 


0.654381 


0.001130 


-0.028649 


0.165992 




0.9 


0.751314 


0.001388 


-0.032590 


0.188004 




1.0 


0.857323 


0.001695 


-0.033125 


0.212903 



We have obtained that the three types of FPs are characterized by the replicon 
eigenvalues Arep coinciding with the corresponding eigenvalues A3 (see tables 1-4 in [5]) 
up to the numerical accuracy. The same identity occurs in the one-loop approximation 
[HI IZ|- Nevertheless, there is no evidence that this identity holds in the higher orders of 
approximation. The obtained values of parameter Arep for the different types of FPs are 
presented in table [TJ The analysis of A (see tables 1-4 of paper [Hj and table [TJ of this 
paper) calculated for the three types of FPs and for different values of p — 1, 2, 3 gives 
the following results. 

(i) 2D and 3D Ising models (d — 2,3; p — 1). Within the one-step RSB subspace 
the only stable FP is the disorder-induced RS FP (type II) with all ReAj > 0, i — 1, 2, 3. 
In the case of the 3D Ising model (d = 3) the disorder-induced RS FP is characterized 
by the complex values Ai and A 2 leading to RG flows which spiral around the FP [TT] 
and giving rise to the oscillating corrections to scaling [Tlj . The type II FP remains 
stable when the continuous RSB modes are taken into account (Arep, Arep > 0). The 
FP of a pure system (type I) and the disorder- induced RSB FP (type III) are unstable 
(A3 ~ Arep < 0). Therefore the critical behavior of weakly disordered 2D and 3D Ising 
systems is realized with the disorder-induced RS FP. 

(ii) 3D XY model (d = 3; p = 2). Although the positive eigenvalues Ai,A2,A3,Afep, 
Arep calculated for the disorder-induced RS FP (type II) indicate that this FP is stable in 
the full RSB space, we have reasons to believe that in the higher orders of approximation 
the type I FP (which corresponds to the critical behavior of a pure system) will become 
stable. Indeed, the small eigenvalues A3, Arep ~ 3 • 10~ 6 indicate the weak stability of the 
disorder- induced RS FP. On the other hand, the two-loop RG analysis of the disordered 
systems without RSB potentials ^2] gives p c = 2.0114 for the borderline between regions 
of stability for the disorder- induced FP (p < p c ) and the FP of a pure system (p > p c ), 
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while the estimation of p c performed in the six-loop approximation using the pseudo 
e expansion gives p c = 1.912 The latter result is in accordance with the Harris 

criterion, which suggests that due to the negative value of the specific heat exponent 
of the pure 3D XY model the critical behaviour of this model should be stable with 
respect to the influence of the quenched disorder at least in the absent of RSB effects. 

(Hi) 3D Heisenberg model (d — 3; p — 3). The disorder- induced RS FP is unstable 
due to the negative eigenvalues A 2 ,A 3 = \i% < 0. The disorder-induced RSB FP is 
stable in the directions corresponding to the continuous RSB modes for 0.1 < xo < 0.2 
(Arep,Arep > 0), but this FP is unstable in the directions corresponding to the one-step 
RSB modes for any xq (A 2 < 0). Additionally, both FPs are characterized by the non- 
physical values of coordinate g* < 0, so that the RS FP of a pure system is the only 
stable FP in the full RSB space. Therefore the influence of weak quenched disorder on 
the critical behaviour of the 3D Heisenberg model is irrelevant even with taking into 
account the RSB effects. 

4. Conclusions 

The RG investigations carried out in the two-loop approximation show the stability of 
the critical behaviour of weakly disordered systems with respect to the introduction 
of potentials with the general Parisi RSB structure. In dilute Ising-like systems, the 
disorder- induced critical behaviour is realized with RS FP. The critical behavior of XY- 
and Heisenberg-like systems is not affected by the weak quenched disorder. 

Recently, the phase diagram of the random temperature Ising ferromagnet was 
studied within the framework of the Gaussian variational approximation [TE] . The 
spin-glass phase separating the paramagnetic and ferromagnetic phases was found. It 
was shown that the transition from paramagnetic to spin-glass state is second order, 
whereas the transition between spin-glass and ferromagnetic states is first order. It 
was also shown that within the considered approximation there is no RSB in the spin- 
glass phase. If it is just the case, we argue that the paramagnetic-spin-glass transition is 
controlled by the traditional RS FP of a weakly disordered system and that slightly above 
this critical point the asymptotic behaviour is characterized by the critical exponents 
obtained previously for the diluted Ising model within the RS RG treatment. However, 
there are some strong arguments against the existence of an intermediate spin-glass 
phase [T7] so that to define the nature of phase transition controlled by the obtained 
stable RS FP further investigation is clearly necessary. 

Acknowledgments 

Support from the Deutsche Forschungsgemeinschaft (SFB 418) is gratefully acknowl- 
edged. I would also like to thank A. Kudlay for a useful discussion. 



Replicon modes and stability of critical behaviour 



10 



References 

[1] Le Doussal P and Giamarchi T 1995 Phys. Rev. Lett. 74 606 

[2] Dotscnko Vik S, Harris A B, Sherrington D and Stinchcombe R B 1995 J. Phys. A: Math. Gen. 
28 3093 

[3] Dotscnko Vik S and Fcldman D E 1995 J. Phys. A: Math. Gen. 28 5183 
[4] Dotscnko Vik S Usp. Fiz. Nauk 1995 165 481 

[5] Prudnikov V V, Prudnikov P V and Fcdorcnko A A 2001 Phys. Rev. B 63 184201 
[6] Wu X T 1998 Physica A 251 309 

[7] De Cesare L and Mercaldo M T 1999 Phys. Rev. B 59 855 

Dc Ccsarc L and Mercaldo M T 1999 Phys. Rev. B 60 2976 
[8] Mezard M, Parisi G and Virasoro M 1997 Spin-Glass Theory and Beyond (Singapore: World 

Scientific) p 476 
[9] Binder K and Young A P 1986 Rev. Mod. Phys. 58 801 
[10] Folk R, Holovatch Yu and Yavors'kii T 2000 Phys. Rev. B 61 15114 
Pakhnin D V and Sokolov A I 2000 Phys. Rev. B 61 15130 
Prudnikov V V, Prudnikov P V and Fedorcnko A A 2000 Phys. Rev. B 62 8777 
Prudnikov V V, Ivanov A V and Fedorenko A A 1997 Sov. Phys. - JETP Lett. 66 835 
[11] Zinn- Justin J 1996 Quantum field theory and critical phenomena (Oxford: Clarendon) p 1008 
[12] Jug G 1983 Phys. Rev. B 27 609 

[13] Holovatch Yu, Blavats'ka V, Dudka M, von Ferber C, Folk R and Yavors'kii T 2002 Int. J. Mod. 

Phys. B 16 4027 
[14] Khmclnitskii D E 1978 Phys. Lett. A 67 59. 

[15] Holovatch Yu, Dudka M and Yavors'kii T 2001 J. Phys. Sudies 5 233 
[16] Tarjus G and Dotsenko V 2002 J. Phys. A: Math. Gen. 35 1627 
[17] Sherrington D 1980 Phys. Rev. B 22 5553 



